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ABSTRACT 

We study the properties of the broad line region in blazars by comparing the virial 
estimate of black hole masses with that derived from the mass of the host galaxies. 
The former is sensitive to the width of broad lines, i.e., to the projection of the velocity 
of line-emitting clouds along the line of sight; the latter is not. This comparison allows 
us to constrain the deprojection factor /, thus revealing general properties of the 
geometry of the broad line region. We show that blazars tend to have 1) higher / values 
than the quasars of our reference sample: (/BLLacs) = 6.9±2.3, (/blazars) = 5.6±1.3and 
(/quasars) = 2.0 ± 0.3; 2) relatively narrow broad emission lines; 3) modest equivalent 
widths, as expected because of the occurrence of jet emission at very low inclination 
angles. In a disc-like sketch of the broad line region, these results indicate a pole-on 
view of a flat geometry in blazars. This consistently extends the orientation-dependent 
unified model of active nuclei to the geometry of the broad line region. 
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1 INTRODUCTION 

The estimate of black hole masses (.Mbh) in distant galaxies 
is feasible only in active galactic nuclei (AGN), where we 
can infer the dynamics of gas close to the black holes (BHs) 
through its bright continuum and line emission. A common 
practice for type-1 AGN (where broad emission lines are 
observed) consists in adopting the virial paradigm: 



The virial technique suffers of a number of limitations. 
In particular, the deprojection value / is not known, and is 
strongly sensitive to the (unknown) geometry and orienta- 
tion of the BLR. In th e commonly assum ed isotropic picture 
of the BLR, (/) = y/z/i llNetzedll990h. Assuming instead 
a disc-like geometry (e.g.. iMcLure fc Dunlop|l2001 ). 



-Mbh = 



-Rblr^Ilr 
G 



/ = 0.5 



(1) 



-1/2 



(3) 



where -Rblr and «blr are respectively the orbital radii and 
the velocity of the gas clouds in the broad line region (BLR), 
and G is the gravitational constant. The cloud velocity is 
usually inferred from the width of broad lines, assuming a 
deprojection factor /: 



«blr = / ■ FWHM 



(2) 



while -Rblr is d erived from the co n tinuum lum i nosity of the 
accre tion disc l|Kaspi et all |2000| . |2005l , 120071 ; iBentz et all 
2009). This allows an est imate of ./VI bh from single-epoch 
spectra: Ishen et all (2010) applied this technique to the huge 
spectroscopic sample of quasars in the Sloan Digital Sky 
Survey (SDSS), thus providing .Mbh for ~ 100, 000 sources. 



where H/R is the aspect ratio of the disc and i? is the angle 
between the normal to the disc and the line of sight. In 
the limit of geometrically thin discs, equation [3] reduces to 
/ = 0.5/(sin$), ranging between 0.5 (edge-on discs) and oo 
(face-on discs). Taking into account the thickness of the disc, 
/ 0.5 (H/R)- 1 when r? -s- For disc-like BLRs (/) « 1 
almost independently of H/R (|Decarli et al.ll2008al h 

From an observational point of view, average values for 
/ have been estimated by comparing the mean offset o f .Mbh 
estimates based on different techniques. lOnken et al.l (|2004l ) 
compared the .Mbh values expected from the .Mbh-o"* rela- 
tion for a set of local type-1 AGN with those obtained from 
the FWHM of KB throu gh equation □ finding (f(KB)) = 
2.fl ICollin et all (2006) repeated the same exercise for a 
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1 The average value reported here has bee n recomputed to ac - 
count for the different definition of / used in lOnken et al.l J2004T) . 
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larger set of objects, and found a dependence of (/(H/3)) on 
the broad line width ({/} « 1.5 for FWHM<2000 km/s, 
(/) w 0.7 for FWHM>6000 km/s), which could be at- 
tribu ted to the Eddingto n factor or to the BLR orienta- 
tion. Decarli et al.l (l2008bl) expanded on the previous results 



bv lLabita et all (2006). comparing virial estimates of A^bh 
based on H/3 and C IV lines with those inferred from the host 
alaxy luminosities assuming the scaling in IBettoni et al.l 
2003h . They found (/(H/3)) w 1.6, (/(Civ)) « 2.4 and a 
strong dependence of / on the widths of broad lines, i.e., 
/ > 2 values are found only in quasars w ith FWHM<5000 
km/s (see also lLa Mura et al.|[2007l . l2009ll . The different val- 
ues of (/) for different broad emission lines can be explained 
by a different geometry of the corre sponding emitting re- 
gio ns (see alsolFine et a l] |2008l , 120101 1 . 

IDecarli et al. I (|2008al ) noticed that the claimed "black 
hole mass defici t" observed in Narrow-L ine Seyfert 1 galax- 
ies (NLSls; see iGrupe fe Mathurl 120041 ) disappears once a 
disc-like BLR is assumed. NLSls show relatively narrow 
permitted lines (500 km/s < FWHM < 2000 km/s), small 
-Mbh with respect to the predictions from the A'fBH-host 
galaxy relations, and high Eddington ratios. Assuming that 
the BLR is fiat and seen almost pole-on, from the statistics 
of NLSls with respect to 'normal' Seyfert 1 galaxies we can 
infer the average / value for NLSls. Interestingly, once cor- 
rected for this higher deprojection values, line widths yield 
to distributions of A^bh and Eddington ratios almost in- 
distinguishable from normal Seyfert 1 galaxies. The strong 
radio loudness wi th blazar-like prop erties observed in a large 
fraction of NLS1 (|Yuan et al.| [2008) and the recent detecti on 
of GeV emission in a handful of NLSls (|Abdo et al.ll2009al lbh 
also reveal the presence of beamed emission almost aligned 
with the line of sight, supporting a pole-on picture of the 
central engine of these AGN. 

Other indications supporting a flat geometry for the 
BLR come from 1) the anti-correlation between the width 
of broad emission lines and the "compactness" parameter 
defined from radio observations as the luminosity ratio be- 
tween the core and the lobe emission in radio galaxies and 
radio-loud quasars dWills fc Brownelll986l ; iBrothertonll 19961 ; 
IVestergaard. Wilkes fc Barthe]||2000l ); 2) the dependence of 
the FWH M of broad lines on the Equivalen t Width of the 
[Om] line (|Risaliti. Salvati fc Marconill2010l '). 

In order to test the fiat geometry picture of the BLR, 
we here extend our previous st udies to blazars. Accord- 
ing to the unified mode l of AGN (|Antonucci fc Millerj[l985l ; 
lUrrv fc Padovamlll995l ). blazar jets (and specifically BL Lac 
objects) are expected to show small inclination angles with 
respect to the line of sight. If the BLR has a flat geometry 
with a rotational axis aligned to the jet, blazars should show 
relatively narrow FWHM of broad lines and high values of 
/. We exploit a newly defined recipe to estimate the ra- 
dius of the BLR from the lum inosity of broad lines (see also 
Vcstcrgaard & Peterson 2006) in order to avoid the contam- 
ination from the jet light in the optical continuum emission. 
Our virial estimates of AIbh are then compared to the es- 
timates we have from the A^BH-A'Ihost relation in order to 
derive /. 

This paper is organized as follows: In section [2] we 
present our sample. Our analysis is described in section 
while in section [4] we present our observational results and 
compare them with expectations from simple BLR models. 



Conclusions are drawn in section [5] Throughout the paper 
we will assume a concordance cosmology with Ho = 70 
km/s/Mpc, fl m — 0.3 and J1a = 0.7. Results from other 
papers are accordingly scaled to this cosmology. 



2 THE SAMPLE 

In order to perform our analysis, we require an estimate 
of the host galaxy luminosity and a characterization of 
the properties (flux, width) of at least one broad line. 
The bla zar sample is most ly derived from the com pila- 
tions by Scarpa etlrtl |2000l ). iKotilainen et al.l l| 19981 ) and 



Nilsson et al 



.20031 ). Few other sources are taken from 
Wurtz et al.P|l996h . iNilsson et al.l i|2009l ) and IDecarli et all 
( 2010bh . We select all the blazars with a reliable estimate 



of the host galaxy luminosity. We adopt published values 
of Mr (host) or convert available bands into R computing 
the fc-correction throu gh the Elliptical galaxy template by 
iMannucci et al.l |200 if ). Out of this selection, we pick up 
only those objects with a broad line detected in an opti- 
cal spectrum. Main data sources for the spectroscopy are 
the zbllac databas^3 and the SDSS archive. For BL Lac 
an d OJ287, we use the FWHM and flu x val ues published 
by Capetti. Raiteri fc Buttiglionel (j2O10T ) and lNilsson et all 
(2010) respectively. It turns out that the only two lines of 
interest are Mg Ha2798 and Ha. The former is the only broad 
line observed at 0.5 < z < 1.6. We preferred the latter rather 
than H/3 because it is ~ 3x brighter. The whole blazar sam- 
ple consists of 18 objects. 

In order to understand the role of orientation of the 
broad line region with respect to 'normal' quasars, we in- 
clude a c omparison sample of q uasars from the data set pre- 
sented in lDecarli et al.l (|2010al lbl). We consider only objects 
where the Mgn or Ha lines have been fully covered in the 
observed wavelength range. The reference sample consists of 
39 sources. Table [T] summarizes the objects in both samples. 



3 DATA ANALYSIS 
3.1 Measured quantities 

With the only exception of BL Lac and OJ287, 
for which published values have been use d 
(|Capetti. Raiteri fc Buttiglionel l20ld : INilsson et all 120091 ') . 
the spectra of all the blazars and reference quasars have 
been fitted with the fitting to o ls deve loped by our group 
and p resented in IDecarli et al . (2010a) and |Pe Rosa et ail 
(2010). Namely, in order to infer the continuum and line 
properties, we apply a multicomponent fit to the continuum 
emission (power-law + Fen model), subtract it, and fit 
the residual with the sum of t wo Gaussian functions with 
the same peak wave length (see IDecarli et al. 2008J], l2010al : 
iDe Rosa et~ai1l2010l). The iron emission is modelled through 
the IVestergaard fc Wilkes! (|200lh template in the UV and 
our own spectrum of I ZwOOl in the optical wavelengths, 
after removing main emission lines and convolving to a 
Gaussian profile to mimic the Doppler broadening. Free 
parameters in the continuum estimate are the slope and 



Sec http: / /www. oapd.inaf.it/zbllac/. 
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Table 1. The sample. (1) Source name. (2—3) Source coordinate in J2000. (4) Redshift. (5) Absolute rest-frame i?-band magnitude of the 
host galaxy. (6) Adop ted Mass to Light ratio, in s olar units. (7) A^bh deri ved from the host ga laxy luminosity. (8) Refe ren ce for the host 
galaxy lu minosity: a- iKotilainen et al.l (l99Sh : b- IScarpa et all j200Ch : c- iNilsson et aP I^OOSt ): d- IWurtz et atl jl996h : e- iDecarli et all 
i2010"tj) : f- lNilsson et alj <l200yt) . 



Name R.A. Dec. z M R (host) M/L log MBH(ima) Ref. 

[mag] [M Q /L ] [M ] 

(1) (2) (3) (4) (5) (6) (7) (8) 



Blazars 
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Table 2. Results from spectroscopy. (1) Source name. (2) Line used to estimate .Mbh- (3) FWHM. (4) Line luminosity. (5) Power-law 
monochromatic luminosity at 3000 A (for Mgn) or 5100 A (for Ha). (6) Line Eq uivalent Width. (7) A ^bh estimate based on the line 
luminosity and width. (8) Geometrical fact or /, (9) Reference f or the spectra: a- [ Decarli et alJ i BoiOal): b- ISbaxufatti et al.l (|2005h : c- 
ISbarufatti et alj<2006ah; d- SPSS arc hive; e-lNilsson et aljfeoid); f-lLabita et al.l <200d) ; e- |Pian. Falomo fc Treves! j2005h : h- HST STIS 
archive; i- ISbarufatti et al.l j2006bl 1; i- ICapetti. Raiteri fe Buttielione! [bold ). 



Name Line FWHM log L linc log \L X E.W. log Mbh (spec) / Ref. 

[km/s] [crg/s] [crg/s] [A] [M ] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



Blazars 



PKS0403-132 


Mgn 


3451 


43.86 


45.77 


47.5 


8.6 


0.6 


a 


PKS0405-123 


Mgn 


3933 


44.49 


46.40 


14.6 


9.0 


0.9 


a 


PKS0420-014 


Mgn 


3607 


43.92 


45.83 


18.8 


8.7 


2.2 


a 


PKS0426-380 


Mgn 


4419 


42.76 


44.67 


3.0 


8.3 


1.7 


b 


PKS0521-36 


Ha 


4657 


41.71 


42.94 


51.9 


7.2 


7.0 


c 


OJ049 


Ha 


2553 


41.55 


42.78 


3.5 


6.6 


21.3 


d 


OJ287 


Ha 


3710 


41.83 


43.06 


0.5 


7.1 


5.6 


c 


3C273 


Mgn 


4074 


43.88 


45.79 


22.4 


8.8 


1.4 


f 


1249+174 


Mgn 


3040 


43.01 


44.91 


27.5 


8.1 


1.0 


d 


3C279 


Mgn 


6534 


43.29 


45.20 


8.6 


8.9 


1.6 


g 


PKS1519-273 


Mgn 


1821 


42.37 


44.28 


0.4 


7.1 


7.9 


b 


4C1460 


Mgn 


2697 


41.86 


43.76 


1.1 


7.5 


9.6 


d 


3C345 


Mgn 


4628 


42.77 


44.68 


13.5 


8.4 


6.0 


a 


3C371 


Ha 


2333 


42.83 


44.07 


14.9 


7.4 


6.5 


h 


PKS2131-021 


Mgn 


3668 


42.61 


44.52 


2.5 


8.1 


1.7 


i 


BL Lac 


Ha 


4600 


41.59 


42.82 


9.9 


7.2 


6.1 


j 


PKS2201+04 


Ha 


3556 


40.61 


41.84 


18.8 


6.3 


15.0 


c 


PKS2345-167 


Mgn 


3068 


43.40 


45.31 


20.3 


8.3 


3.9 


a 



Reference quasar sample 



PKS0000-177 


Mgn 


6475 


43.74 


45.65 


19.9 


9.1 


1.0 


a 


SGP5:46 


Mgn 


3504 


42.94 


44.85 


23.4 


8.2 


1.4 


a 


0054+144 


Ha 


8635 


42.53 


43.76 


399.0 


8.3 


2.4 


a 


PKS0100-270 


Mgn 


3384 


44.22 


46.12 


15.5 


8.8 


0.9 


a 


0119-370 


Mgn 


6229 


43.58 


45.49 


60.3 


9.0 


0.9 


a 


0133+207 


Mgn 


12408 


43.72 


45.62 


77.4 


9.7 


0.4 


a 


0152-4055 


Mgn 


4924 


43.91 


45.82 


31.1 


9.0 


0.7 


a 


PKS0155-495 


Mgn 


5395 


43.81 


45.72 


50.9 


9.0 


1.2 


a 


PKS0159-11 


Mgn 


4140 


44.17 


46.08 


26.4 


8.9 


0.7 


a 


PB6708 


Mgn 


3156 


43.58 


45.48 


53.8 


8.4 


1.4 


a 


KUV03086-0447 


Mgn 


2441 


43.94 


45.84 


19.0 


8.4 


2.1 


a 


US3828 


Mgn 


3779 


43.95 


45.85 


47.2 


8.7 


1.1 


a 


PKS0348-120 


Mgn 


3699 


43.92 


45.83 


27.3 


8.7 


1.8 


a 


PKS0349-14 


Mgn 


5841 


44.21 


46.12 


21.6 


9.2 


1.9 


a 


PKS0402-362 


Mgn 


4619 


44.19 


46.10 


17.8 


9.0 


1.3 


a 


PKS0414-06 


Mgn 


5732 


44.38 


46.29 


45.7 


9.3 


1.5 


a 


PKS0440-00 


Mgn 


5585 


43.29 


45.20 


23.7 


8.8 


1.3 


a 


0624+6907 


Mgn 


3989 


44.50 


46.40 


24.5 


9.0 


2.0 


a 


PKS0710+11 


Mgn 


11759 


44.30 


46.21 


37.7 


9.9 


1.0 


a 


MS07546+3928 


Ha 


2126 


43.52 


44.75 


362.7 


7.8 


3.1 


a 


MS0824.2+0327 


Mgn 


2812 


43.76 


45.66 


29.6 


8.4 


1.7 


a 


MS08287+6614 


Mgn 


6403 


43.54 


45.44 


32.8 


9.0 


0.9 


a 


PKS0838+13 


Mgn 


4383 


43.42 


45.32 


132.1 


8.6 


1.4 


a 


US1867 


Mgn 


2693 


43.64 


45.55 


27.2 


8.3 


7.7 


a 


TON392 


Mgn 


4899 


44.34 


46.25 


31.5 


9.2 


0.9 


a 


MS09441+1333 


Ha 


2564 


43.12 


44.35 


344.6 


7.7 


6.2 


a 


1004+130 


Ha 


6356 


43.72 


44.95 


238.9 


8.8 


1.7 


a 


Z101733-0203 


Mgn 


2914 


42.61 


44.52 


28.5 


7.9 


1.3 


a 


PKS1015-31 


Mgn 


3835 


43.14 


45.04 


1001.3 


8.4 


2.5 


a 


PKS1018-42 


Mgn 


4105 


44.18 


46.09 


34.8 


8.9 


1.8 


a 


1150+497 


Ha 


3523 


43.47 


44.70 


439.4 


8.2 


3.2 


a 


1402+436 


Ha 


2556 


43.01 


44.24 


249.5 


7.6 


6.1 


a 


1444+407 


Ha 


2694 


42.69 


43.92 


355.2 


7.1 


4.8 


a 


PKSJ1511-10 


Mgn 


3285 


43.72 


45.63 


16.0 


8.5 


1.3 


a 


HS1623+7313 


Mgn 


3250 


43.38 


45.29 


36.3 


8.3 


1.3 


a 


MC2112+172 


Mgn 


4221 


43.28 


45.19 


38.2 


8.5 


2.3 


a 


Q2225-403B 


Mgn 


4479 


43.22 


45.13 


48.2 


8.5 


1.4 


a 


2247+140 


Ha 


2832 


42.86 


44.10 


120.8 


7.6 


5.6 


a 
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Figure 1. Examples of the continuum and line fits. The spectra of 
two blazars (PKS0426-380 and OJ049) and two quasars (US3828 
and 1004+130) are shown. The solid black lines show the observed 
spectra, the solid red lines mark the model, which results from the 
superposition of a power-law continuum (dotted black lines), the 
iron and host galaxy templates (dotted light blue lines) and a line 
model (dotted red lines). The residuals after model subtraction 
are shown as blue dashed lines (a colour version of this plot is 
available in the electronic edition of the journal). 



normalization of the power-law and the intensity of the Fe II 
feature. The iron multiplets are found to be negligible in the 
wavelength range bracketing Ha. Two sources (PKS0521-36 
and PKS2201+04) show significant emission from the host 
galaxy in their optical spectra. We modelled this componen t 
using the Elliptical template by iMannucci et al.l i|200 ll ). 
According to our best fit, at 5100 A, the host galaxies 
of these two objects account for ~ 30% of the specific 
flux. Two more objects (3C371 and OJ049) also show 
marginal (< 5 % of the flux density at 5100 A; see Figure 
[TJ contribution from the host galaxies. 

The line fit has 5 free parameters, namely the line peak, 
and the width and amplitude of the two Gaussian curves. 
Figure [T] shows four examples of fitting. From the line model 
we compute the Full Width at Half Maximum (FWHM), the 
integrated line fl ux (fline), and the li ne Equivalent Width 
(EW). Following |Pe Rosa et all l|2010] ). we compute uncer- 
tainties in the parameters by marginalizing the distribution 
of values yielding \ 2 < Xbest + 1- Fitting errors in the lumi- 
nosities (both for the line and the continuum) are usually 
~ 5%, and never exceed 10%. Typ ical errors on the FW HM 
are ~ 10% (see also the analysis in lDecarli et al.|[2008bl ). We 
remark that these errors do not contribute significantly to 
the uncertainties on AIbh and /, which are dominated by 
the scatter in the i?BLR-luminosity relations (~ 0.4 dex). 



3.2 Derived quantities 

The J?BLR-himinosity relations quoted above refer to the 
continuum emission of the accretion disc. This quantity is 
not directly accessible in blazars, where significant contam- 
ination from the synchrotron emission of the jet may occur. 
However, quasar spectra show a strong correlation between 
the luminosity of broad lines and the continuum luminosity 
of the accretion disc. This is not unexpected, as broad lines 
are emitted by gas direc t ly ph otoionized by the continuum. 
IVestergaard fe Peterson! (|2006l ) used the broad H/3 line lu- 
minosity as a proxy for the BLR radius in quasars. In the 
present case, this approach benefits of avoiding contamina- 
tions by the jet emission, hence representing a viable tool to 
estimate -Rblr even in blazars. 

In Figure [2] we plot the c ontinuum v e rsus t he broad line 
luminosities in quasars from IShen et al l (|2010l ) for Ha and 
Mgn. The latter shows a clear 1-to-l correlation spanning 
several order of magnitudes. The average ratio between the 
continuum and the line luminosities is: 



log 



AL A (3000A) 
L linc (MgII) 



= 1.91 ±0.26 



(4) 



where the error is the rms, and it was computed over 85515 
quasars. Concerning Ha, the scene is more complex. At low 
luminosities (< 10 44 erg/s in the continuum) , significant 
contam ination from the host galaxy is present. iBentz et al.l 
(2009) showed that the removal of this contribution signifi- 
cantly improves the quality of the J?BLR-luminosity relation 
in the (rest frame) optical wavelengths. We therefore focus 
only on the bright end of the continuum luminosity distribu- 
tion, where the role of starlight is negligible. A similar cut is 
also applied to the line luminosity, since fits of faint lines are 
less reliable due to the contaminations from narrow Ha and 
[Nil] lines. The average and rms values of the continuum to 
line luminosity ratios (computed over 1162 objects) is: 



log 



AL A (5100,4) 
Li inc (Ha) 



1.23 ±0.14. 



(5) 



We used equations [4] and [5] to infer the disc contribution to 
the continuum luminosity. Since only one of the two lines 
can be observed in optical spectra (Mgn or Ha, according 
to the target redshift), we adopted equation [4] for objects at 
z > 0.5 and equation [5] for the remaining sources. 

Once the line width, luminosity and the (corrected) con- 
tinuum lum inosity are known , we u se th e J^BLR-lum i nosity 
relations bv lMcLure fe Jarvisl (|2002l ) and lKaspi et all <|2007h 
to compute Mbh through equations Q~H5] These estimates 
are compared with the A^bh values we derive from the host 
galaxy. We convert host galaxy luminosities into masses by 
assuming that the stellar populations are old, generated in 
a singl e burst of instantane ous star formation at high-z (fol- 
lowing Dccarli et al . 2010b, we assume Zburst = 5), then pas- 
sively evolving. We adopt the mass-to -light ratio evolution 
predicted bv lBruzual fe Chariot! |2003l ) . The black hole mass 



3 All the sources i n our study have e lliptical host galaxie s 
iScarpa et al.l l200d : iNilsson et ail 120031 : iDecarli et ail fooiOaT) . 
Therefore, we do not discriminate between the bulge and the 
whole galaxy. This is also consistent with the assumption of an 
old stellar population. 
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Figure 2. The continuum to line luminosity ratios for Ha (left panel) and Mgn (right panel). Data are taken from' Shen et al 
For the sake of clarity, in the Mgn plots, only 1 object out of 5 is shown. Contours enclosing > 10, 20, 40, 80, . . .objects are also plotted. 
Dash lines show the average ratios. The plot region used to compute it for Ha is marked with dotted lines. The faint end of the continuum 
luminosity distribution in the optical wavelengths is clearly deviating from the 1-to-l relation, because of the contribution of the host 
galaxy starlight. The same is not true in the UV range, where the host galaxy contribution is negligible throughout the luminosity range 
of interest, and a clear linear correlation appears. 



is then e stimated through the (redshift-dependent) relation 
found in lDecarli et al.1 (|2010bl ): 

Mp,H 

log = (0.28 ± 0.06)z - (2.91 ± 0.06) (6) 

bulge 

We therefore can estimate / from the ratio between the 
image-based estimate of Mbh and the one obtained from 
equation [T] using FWHM instead of ublr: 

_ / Mbu (image) _ / GMbh (image) 

V XBH(virial) ~ V -RblrFWHM 2 ( ' 

Tables[T]and[2]list the main measured and derived quan- 
tities for the blazar and reference samples. 



4 RESULTS AND DISCUSSION 

In Figure[3]we compare the Mbh estimates obtained in the 
virial assumption, assuming / = 1, with the bulge mass in- 
ferred from the host galaxy luminosity. It is apparent that 
blazars, and specifically BL Lac objects, show smaller virial 
A4bh values than quasars for a given host galaxy mass, i.e., 
higher values of / are required for blazars. We note that the 
redshift-dependent evolution of the AlBH/.Mhost scaling re- 
lation has not been applied at this stage. The bulge mass 
range of blazars and quasars sampled in this study is the 
same, so no explicit bias due to the host galaxy is expected 
to affect our results. The same is true if we divide our sample 
according to the broad line width: No obvious segregation 
in terms of bulge mass is observed. Most importantly, when 



we pick up pure BL Lac objects (with broad line Equiv- 
alent Width EW<10 we automatically select blazars 
among those showing the largest deviations between A4bh 
and A^buigc, as expected if we assume that BL Lac objects 
show the smallest inclination angles amongst blazars. 

Figure U shows the / v alues as a function o f the line 
FWHM. As already noted in lDecarli et al.1 (|2008bh , an anti- 
correlation between / and FWHM is apparent. This sup- 
ports the idea of a flat geometry viewed at different inclina- 
tion angles. In particular, only few objects with / > 2 have 
FWHM>4000 km/s, and no one of them has FWHM>5000 
km/s. It is also apparent that blazars have, on average, 
higher / than quasars. Namely, (/BLLacs) = 6.9 ± 2.3, 
(/blazars) = 5.6 ± 1.3 and (/quasars) = 2.0 ± 0.3, with all 
the BL Lac objects having / > 2. All these indications sup- 
port a flat geometry of the BLR, with smaller inclination 
angles for blazars and BL Lacs in particular. 

Further evidence of this appears from the distribution 
of sources in the /-EW plane. As shown in Figure [2] the 
luminosity ratio between broad lines and continua is almost 
consta nt in the average quasar populatio n (see also the anal- 
ysis in lRisaliti. Salvati fc Marconil l2010). In BL Lac objects, 
the continuum emission from the jet contaminates the op- 
tical spectrum, so that the ratio between the line flux and 



4 Here we slightly relaxed the usual EW<5 A threshold in order 
to include BL Lac and 3C279. We note that, in order to study 
broad emission lines, we are biased towards spectra collected dur- 
ing relatively faint states of the continuum emission of BL Lac 
objects. 
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Figure 3. Comparison between black hole masses (derived 
through the virial technique, assuming FWHM=i>blr) an d the 
mass of the host galaxy bulge. The local iMarconi &: Hunt! ||2003T ) 
relation is plotted as dashed lines. Different symbols mark differ- 
ent subsamples: Filled / empty circles mark the blazars / refer- 
ence quasars. Filled / empty squares mark objects with line width 
smaller / larger than 4000 km/s. Filled / empty triangles mark 
objects with EW below / exceeding 10 A. Blazars (and BL Lac 
objects in particular) show smaller virial jMbh than expected 
from the host galaxy mass, i.e., larger / values. 



Figure 4. The geometrical factor / as a function of the width of 
the lines (FWHM). Dotted lines show the expected values of / 
for a thin disc model of the BLR, with different inclination angles 
(labelled in degrees). Diagonal, dashed lines mark the loci at fixed 
cloud velocities in the thin disc model, namely t)BLR=5000, 10000 
and 20000 km/s. Symbol code is the same as in Figure [3] Blazars 
(and specifically BL Lac objects) show relatively narrow lines and 
high / values, consistently with the picture of a flat BLR geometry 
observed pole-on. 



the measured continuum (i.e., the line EW) is decreased. In 
Figure[5]we show that a) no quasar is found with small EW, 
as expected; b) the higher is /, the smaller is the EW. In 
particular, all the objects with EW<10 A have / > 2. This 
picture is in good agreement with the expectations from the 
unified model of AGN if we assume a flat geometry for the 
BLR. 

In order to put these results in the context of a quan- 
titatively defined picture of the BLR, we compare the dis- 
tribution of / observed in our study with the expectations 
from equation [3] through Montecarlo realizations. We gen- 
erate a compilation of 10 5 objects with isotropically dis- 
tributed inclination angles i? between the BLR rotational 
axis and the line of sight. We consider as pole-on systems 
those having $ < $ cut = 10°. The cutoff angle is set by 
the typical opening angle of a jet in BL L ac objects (see 
Ghiscllini ct al . 1993; U rry fc Padovani|[l995l ). i.e., we select 
sources in which we expect the jet to contaminate the con- 
tinuum emission in the optical wavelengths. Objects with 
#cut < t? < $max ~ 40° are considered normal type-1 AGN. 
The cutoff between these two classes is smoothed by adding 
an error ay — 10° to the intrinsic values of This un- 
certainty accounts for errors in the estimate of /, intrin- 
sic differences in the opening angles of the jets, and physi- 
cal mis-alignment between the jets and the BLR rotational 
axis. This last effect is expect ed in a disc warped by, e.g ., 
the Bardeen-Petterson effect IjBardeen fc Petterson|[l975l ). 



In particular, for very massive BHs as those discussed here, 
the warping radius of a disc-like distri bution of gas i s com - 
parable with the radius of the BLR (Pe rego et al.l f2009). 
In this picture, an average misalignment of ~ 10° is ex- 
pect ed if the BH is fu elled by a larger scale circumnuclear 
disc (|Dotti et al.ll201ol ). Other physical effects not considered 
her e for sake of simplicity (e.g. self-irradiation driven warps, 
see lPringle 1996; M alonev. Begelman fc Pringldll996l ) could 
also play an additional role in smoothing the cutoff. 

For each source we compute / through equation [3] and 
we compare the resulting distributions with what observed 
for objects with broad line Equivalent Width below 10 A. 
Expected and observed distributions of / are shown in Fig- 
ure [S] where we assumed H/R = 0. Despite the modest 
number of BL Lac objects, the agreement between the pre- 
dicted and observed distributions is remarkably good. The 
Kolmogorov-Smirnov test indicates that the probability that 
the observed distribution of / for objects with EW<10 A is 
obtained by the same parent distribution that generated the 
pole-on sample in the Montecarlo realization is high (40%). 

Concerning H/R, the observation of few, very-high val- 
ues of / (> 10) suggests that, at least for these objects, 
H/R < 0.05. Re-running our Montecarlo realizations with 
H/R = 0.05 yields consistent results: In particular, the prob- 
ability that the pole-on sample and the observed BL Lac 
objects have the same parent distributions is 20%. 
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Figure 5. The geometrical factor / as a function of the Equiva- 
lent Width of broad lines. Symbol code is the same as in Figure 
\3\ A trend is apparent, so that objects with strong jet contamina- 
tion (= small line Equivalent Width) all present large geometrical 
factor values (consistent with a pole-on view of a flat geometry). 




Geometrical factor 



Figure 6. Observed distributions of the geometrical factor / for 
objects with broad line Equivalent Widths exceeding / below 10 A 
(solid empty / shaded histograms), as compared with the output 
of our Montecarlo realizations for pole-on sources (dotted shaded 
histogram; see the text for details). The agreement is remarkable, 
given the modest statistics of BL Lac objects in our sample. 



5 CONCLUSIONS 

We computed virial estimates of A^bh in a set of 18 blazars 
and 39 reference quasars, and compared them with the AIbh 
estimates derived from the host galaxy luminosities. Virial 
black hole masses are derived through broad line luminosi- 
ties and widths, using a newly defined relation between the 
line and the continuum luminosity in quasars. From the com- 
parison of the two Mbh estimates we infer the geometrical 
factor /, responsible for the deprojection of the virial ve- 
locity from the line of sight, and we used / to constrain the 
properties of the broad line region in blazars. We found that: 

i- Blazars, and BL Lac objects in particular, have higher 
/ values on average with respect to the reference quasar 
sample, and both these classes have average / values well 
exceeding the predictions for an isotropic BLR. The average 
values of / for each class of objects are: (/BLLacs) = 6.9±2.3, 

(/blazars) = 5.6 ± 1.3 and (/quasars) = 2.0 ± 0.3. 

ii- The geometrical factor clearly depends on the width 
of broad lines, in the sense that only few objects with 
/ > 2 have FWHM>4000 km/s, and no one of them has 
FWHM>5000 km/s. All the blazars in our sample have rel- 
atively narrow lines (<5000 km/s). 

Hi- An anti-correlation between the geometrical factor and 
the broad line equivalent width is also found, in the sense 
that all the objects with EW<10 A have geometrical factors 
/>2. 

These results rule out an isotropic sketch of the broad line re- 
gion. The observed distribution of / is in excellent agreement 
with the expectations for a geometrically-thin (H/R < 0.05) 
disc picture. In this scheme, the broad line region of BL Lac 
objects is seen almost pole-on. Assuming small or negligible 
tilts (< few degrees) between the disc rotational axis and 
the AGN jet, this conclusion is in full agreement with the 
Unified Model of AGN. 
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